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Abstract 

We calculate exactly the first cumulants of the free energy of a directed polymer in a random medium for the 
geometry of a cylinder. By using the fact that the n-th moment (Z n ) of the partition function is given by 
the ground state energy of a quantum problem of n interacting particles on a ring of length L, we write an 
integral equation allowing to expand these moments in powers of the strength of the disorder 7 or in powers 
of n. For n small and n ~ (Z^) -1 / 2 , the moments (Z n ) take a scaling form which allows to describe all 
the fluctuations of order 1/L of the free energy per unit length of the directed polymer. The distribution of 
these fluctuations is the same as the one found recently in the asymmetric exclusion process, indicating that 
it is characteristic of all the systems described by the Kardar-Parisi-Zhang equation in 1 + 1 dimensions. 

P.A.C.S. numbers: 05.30, 05.70, 64.60 Cn. 

1 Introduct IOH an interesting system from the point of view of non- 
equilibrium phenomena: through the Kardar-Parisi- 

Directed polymers in a random medium is one of the Zhang (KPZ) equation^ it is related to ballistic 

simplest systems for which the effect of strong disor- growth models and, in 1 + 1 dimensions, to the asym- 

der can be studied]!], § ||. At the mean field level, metric simple exclusion process (ASEP)|j, |9). 

it possesses a low temperature phase, with a bro- In the theory of disordered systems, the replica ap- 

ken symmetry of replica^, |f| similar to mean field proach plays a very special role. On the one hand, it 

spin glasses Q. The problem is however much better is one of the most powerful theoretical tools and often 

understood than spin glasses; in particular one can the only possible approach to study some strongly 

write|^, H closed expressions of the mean field free disordered systems. On the other hand it is diffi- 

energy and one can predict the existence of phase cult to tell in advance whether the predictions of the 

transitions in all dimensions d + 1 > 2 + 1. It is also replica approach are correct or not. When it does not 
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work, one can always try to break the symmetry of 
replica^): this usually makes the calculations much 
more complicated without being certain that the re- 
sults become correct. In the replica approach, the 
calculation usually starts with an integer number n 
of replica. Then, as the limit of physical interest is 
the limit n — > 0, one has to extend to non-integer n 
results obtained for integer n. This is in fact the big 
difficulty of the replica approach, so it is useful to 
look at simple examples for which the n dependence 
can be studied in detail. 

This is one of the motivations of the present work, 
where we show how to calculate integer and non- 
integer moments (Z n ) of the partition function Z of a 
directed polymer in 1 + 1 dimensions. The geometry 
we consider is a cylinder infinite in the t direction and 
periodic, of size L, in the x direction (i.e. x + L = x). 
The partition function Z(x, t) of a directed polymer 
joining the points (0,0) and (x,t) on this cylinder is 
given by the path integral 



Z{x,t) 



(x.t) 

T>y(s) exp 

(0,0) 



1 f dy(s) 

2 I ds 



(i) 



where the random medium is characterised by a 
Gaussian white noise rj(x, t) 

( V (x,t)r)(x',t')}= 7 5(x-x') 5(t-t'). (2) 

One of the main goals of the present work is to cal- 
culate the cumulants lim t ^ 00 (ln fe Z(t)) c /t of the free 
energy per unit length of the directed polymer. These 
cumulants are the coefficients of the small n expan- 
sion of E(n, L, 7) defined as 



E(n,L,j) = — lim — In 

t^oo t 



(Z n (x,t)) 



(Z{x,t))> 



(3) 



This E(n, L, 7) was calculated exactly by KardarJlTJ 
for integer n and L = 00. His closed expres- 
sion E(n, 00,7) = —n(n 2 — l)7 2 /24 cannot how- 
ever be continued to all values of n, in particu- 
lar to negative n, as it would violate the fact that 
d 2 E(n, L, 7)/ dn 2 is negative. Therefore one does not 
know the range of validity of this expression. 



The second motivation of the present work is to test 
the universality class of the KPZ equation. The prob- 
lem (|l|) of a directed polymer in a random medium 
is described by the KPZ equation as several other 
problems such as growing interfaces or exclusion 
processes H. For certain models of this class, the 
asymmetric exclusion processes, the distribution of 
the total current Y t integrated over time t has been 
calculated exactly Jll], [T^, [TJ|, [l4[ |l5) in the long time 
limit. For large t, the generating function of this 
integrated current Y t on a ring of L sites takes the 
formfH @ 



ln(e Qyt ) 



A ma x(a:)t, 



(4) 



and it was shown [p. l[ |12|, |13|, [L4j, when L is large 
and when the parameter a in (||) is of order L -3 / 2 , 
that A max (a) takes the following scaling form 



A max (a) - aKx = K 2 G(aK 3 ) 



(5) 



where K\, Ki and if 3 are three constants which de- 
pend on the system size L, the density of particles 
and the asymmetry. 

The interesting aspect of (|^) is that the function 
G(/3) is universal|TJ, Q |TJ| in the sense that it does 
not depend on any of the microscopic parameters 
which define the model. It is given (in a paramet- 
ric form) by 



G((3) 



E 

+00 



p 3/2- 



E ^5/2 ' 



(6) 
(7) 



In the correspondence H between the directed 
polymer problem and the asymmetric exclusion pro- 
cess through the KPZ equation, the role played by 
ln(Z(i)) is the ratio Y t /L. Comparing (exp(ay t )) 
and (Z n (t)) in equations (^, ^), we see that n cor- 
responds to aL and E(n,L,"f) to A max (a). If the 
function G(j3) is characteristic of systems described 
by the KPZ equation, we expect in the scaling regime 
(large L and n ~ L -1 / 2 ), a relation similar to (||) 
between E(n,L,j) (defined by (|^)) and n. This is 
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indeed one of the main results of the present work: 
when L is large and n ~ L -1 / 2 , we find 



E(n, L, 7) = — — 



G(-n^2nIsy). (8) 



It is clear that in order to establish this relation we 
have to calculate non-integer moments of the parti- 
tion function. 

The paper is organised as follows. In section g, 
we recall how the replica approach of ([!]) can be for- 
mulated as a quantum problem with n particles on 
a ring and how this problem can be solved by the 
Bethe ansatz when the noise is 8 correlated as in (g) . 
In section |3|, we write an integral equation (26 ) which, 



together with some symmetry conditions (|2J 
lows to solve the Bethe equations of section £2[ 
main advantage of (|2^) is that the strength c of the 
disorder (where c = jL/2) and the number of replica 
appear as continuous parameters. We show how ex- 
pansions in powers of c or in powers of the number n 
of replica can be obtained from this integral equation. 
In the expansion of the energy E(n, L, 7) in powers 
of c, all the coefficients are polynomials in n. This 
allows us to define E(n,L,j) for a non-integer n at 
least perturbatively in c. At the end of section 
we show how to generate a small n expansion which 
solves the integral equation (^6|). We also give ex- 
plicit expressions up to order n 3 and we notice that 
in this small n expansion of the energy, we have to 
deal with coefficients that are functions of c with a 
zero radius of convergence. The content of sections || 
and H is essentially a recall of a method developed in 
our previous workJl7|. In section |], we show that the 
recursion of section^ which generates all the terms of 
the small n expansion simplifies greatly in the scaling 
regime (c large and n ~ c" 1 / 2 ) allowing to calculate 
all the terms of the expansion and to establish (ft). 



2 A quantum system of n par- 
ticles with £ interactions 

Let us start with a case slightly more general than (^|) 
where the noise rj(x,t) in (Q) is a Gaussian noise 8- 
correlated in time but with some given correlation v 



m space 

(r){x, t)r](x',t')) =7»(i- x') 8(t - t'). (9) 

If we consider the correlation function 
(Z(x\, t) Z(x2, t) ...Z(x n ,t)) of the partition 
function Z(x,t) at points x\, x%, ... , x n , one can 
check H from ||) that it satisfies 

■^-(Z(xi,t)Z(x 2 ,t) . . . Z(x n ,t)) = 
at 

- n (z(x 1 ,t)z(x 2 ,t)...z(x n ,t)). (10) 

where the Hamiltonian TL is given by 

(ii) 



a</3 



and where, because of the cylinder geometry in the 
directed polymer problem, we have x a = x a + L for 
1 < a < n. 

This implies that in the long time limit, 

(Z(x u t)Z(x 2 ,t)...Z(x n ,t))~e- t ^ n > L <"<\ (12) 

where E(n,L,~f) is the ground state energy of ([Tl|) . 

If one takes the limit v(x — x') — > 8(x — x'), 
the energy E(n,L,j) becomes infinite because of 
the constant part nv(0)/2 in (pi]). This diver- 
gence disappears, however, if we consider the ratio 
(Z(x 1 ,t)Z(x 2 ,t)...Z(x n ,t))/Yl a (Z(x a7 t)), and one 
can see that in the long time limit, 

(Z(x 1 ,t)Z(x 2 ,t)...Z(x n ,t)) -tE(n,L n ) (n] 

(Z(x u t)){Z(x 2 ,t))...(Z(x n ,t)) ~ ' 1 1 

where E(n, L, 7) is the ground state energy of the 
Hamiltonian 



n 



^ dxl 



7E 

a<(3 



8(x a - Xp), 



(14) 



where the positions x a of the n particles are on a ring 
of length L. 

Lieb and Liniger have shown that the Bethe 
ansatz allows to calculate the ground state energy 
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E(n, L, 7) of this one dimensional quantum Hamilto- 
nian exactly || £H], |§ || The Bethe 

ansatz consists in looking for a ground state wave 
function 9(xt, . . . ,x n ) of (|l4"|) of the form 



. . . , x n ) = qp e 



2(gi£c P(1) H f-9„Kp(„))/i 



in the region < x\ < • • • < x n < L. The sum in ( |l5| ) 
runs over all the permutations P of {1, . . . , n} and the 
value of ^ in other regions can be deduced from (15) 
by symmetries. One can show]22], |23| , |2^, [l7| ] that (15) 
is the ground state wave function of (|14|) at energy 



E(n,L n ) = -- ^ 

1<Q<?7 



(16) 



if the q a are the solutions of the n coupled equations 
^2g a _ TT la — Qp +c 



n 



VP 



(17) 



(18) 



obtained by continuity from the solution {q a } = {0} 
at c — where 

jL 

Moreover, the q a are all different and the ground state 
is symmetric ({q a } = {—Qa})- (See for instance p^| . 
Note that ikj and c in are here j-g, and —7; so 
our c defined by (|l^) and the c in f22|| are different.) 
If we introduce the polynomial P(^0 



P(X) = l[(X-q a ), 

the system of equations ( |T7j ) becomes 
e^Pfe - c) + e- q "P(q a + c) 



0, 



(19) 



(20) 



for any 1 < a < n, and we have from the symmetry 
of the ground state 



P(-X) = (-1)"P(X). 



(21) 



The knowledge of the polynomial P(X) determines 
the energy (|l6|) as 



P(X) = X n 




(22) 



(using (|T^) and the fact that la = 0.) 

For small c, it is possible to solve directly (|2(]) and 
to determine the q a (see appendix ^). This leads 
to the following expression of the ground state en- 
ergy (|l|) 



( 15 ) E(n,L,j) 



180 



(23) 



We see that the first coefficients of the small c ex- 
pansion are polynomial in n. In fact, following the 
approach of appendix [5], one can see that each coef- 
ficient of the small c expansion of E{n, L, 7) is poly- 
nomial in n, allowing to define, at least perturba- 
tively in c, the ground state energy E(n, L, 7) for 
non-integer n. The approach of appendix |D| becomes 
however quickly complicated. This is why in the next 
section we develop a different approach jl7| based on 
the integral equation (Uw. 



3 Solution of the Bethe ansatz 
using an integral equation 

In this section we recall the approach developed in 
our previous work]!?]], which consists in writing an 
integral equation where c and n appear as continuous 
parameters and which allows to expand the energy in 
powers of c as well as in powers of n. 

Let us introduce the following function of {q a }'. 



B{ u ) = -e^^^Y p(q a )e^ u -^ 



where the parameters p{q a ) are defined by 



(24) 



p(Qa) = n 



gq - qp + C 

q a -qp 



(25) 



If the {q a } are given by the solution of ( |l7j ) which 
corresponds to the ground state, one can show (see 
appendix |A|) that the function B(u) satisfies the in- 
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tcgral equation 

B{1 +u)- B(l - it) = 

nc f dv e -c(v 2 -uv)/2 B ^ _ v ^b(\ + u -v). 
Jo 

and the following two conditions 

5(1) = 1, 
B{u) = B(-u). 



(26) 



(27) 
(28) 



Moreover, the energy jiff ) can be extracted from the 
knowledge of B(u) through 



12 2 V ' 



(29) 



The derivation of |27], |2g, g9|) is given in ap- 
pendix |A|. We are now going to see how one can find 



pcrturbatively in c or in n the solution of ( 26 
and, consequently, the ground state energy 



27], 



for any u. We put this back into (|26|) and we get at 
first order in c 



Bi(l + u) - Bi(l - u) = ™. 



(33) 



Again, there is a unique polynomial solution which 
satisfies the facts that Bi(u) is even and that B\(l) = 
0: 



Bl (u) = -(u 2 -l). 



(34) 



It is easy to see from (26) that at any order in c, we 
have to solve 



B k (l + u) - B k (l - u) =<j>k(u), 



(35) 



where </>fc(it) is a polynomial odd in u. There is a 
unique even polynomial Bk(u) solution of ( |35| ) satis- 
fying -Bfc(l) = 0: it is one degree higher than 4>k(u) 
and can be determined by equating each power of u 
on both sides of (|35|) . [ Alternatively, we found a way 
of writing the solution for any <pk(u): 



3.1 Expansion in powers of c 

To obtain the small c expansion of B(u) for arbi- 
trary n, we write 



B(u) = B (u) + cBx{u) + c 2 B 2 (u) + . . 



(30) 



Conditions (|27j) and (|2J) impose that B (0) = 1 and 
all -Bfc(l) = for k > 0, and that the B/,(u) are all 
even. Moreover, as can be seen directly from ( p7|), 
the q a scale like yfc when c is small. (Appendix p| 
shows how to obtain the small c expansion of the q a .) 
This implies from the definition ( p4] ) of B(u) that all 
the Bk{u) are polynomials in u. 

At zcro-th order in c, (Eq) becomes: 



B (l + u) -B (l-u) =0. 



(31) 



The only polynomial solution of ( pl|) consistent 
with (H, ||), i.e. B {u) = B {-u) and B (l) = 1 is 
simply 



B (u) 



(32) 



B k (u) 



f Jdv4> k {v) + S -±{<t>' k {u)-<t>'k{-L)) 



+1 (^- i) M-^r i) d)) + 



/2 
(36) 



where the Sfc are the coefficients of the expansion of 
a;/ sinhx in powers of x (i.e. as xj sinhx = 1 — x 2 /6 + 
7x 4 /360+. . . , one has sq — 1, si — —1/6, s 2 = 7/360, 
...).] 

This procedure gives for the first terms 

cn(u 2 - 1) c 2 n(2n + l)(u 2 - l) 2 



B(u) = 1 



96 

4n(: 
+2(u 2 - 3) 



. ., . , ,'5n 2 (u 2 -l) + 4n(2u 2 - 1) 



5760 



+ 0(c 4 ). 



(37) 
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The energy can then be deduced from 



E{n,L,i) 



, n(n-l) 
J L 2 



c 



n 



2 12 180 

2 



(38) 



n 



n 



1512 1260 



c 4 + .. 



(For (|38j), we used more terms than given above in 
B(u).) Of course, this expression agrees with ( p3f ) 
obtained directly by expanding the q a . 

3.2 Expansion in powers of n 

The number of particles n is a priori an integer. How- 
ever, when we look at the small c expansion ( |37| ) of 
B(u) or ( |38j ) of the energy, we see that at any given 
order in c the expression is polynomial in n. There- 
fore, one can extend the definition of the small c ex- 
pansion of B(u) or of E(n, L, 7) to non-integer n. We 
can also collect in the small c expansion of B(u) all 
the terms proportional to n and call this series b\ (u). 
From @ we see that 



bi (u) 



(w 2 - 1) (u 2 - l) 2 



4 96 
(m 2 -1) 2 (m 2 -3) 
2880 



0(c 4 



(39) 



More generally, we can collect all the terms propor- 
tional to n k in the small c expansion and call the 
series bk{u). This means that we can write B(u) as 
a power series in n 



B(u) = 1 + nbi(u) + n 2 b 2 (u) + 



(40) 



where all the b k (u) are defined perturbatively in c. 
Conditions ( p7[ p8| ) impose that all the bk(u) are even 
and that b k (l) = for all fc > 1. We define b a (u) = 
1 for consistency. (It is easy to see in the small c 
expansion that if n = 0, then -B(it) = 1.) 

We are now going to describe the procedure we 
used|l7j] to determine the whole function b\{u) and 
eventually all the bk(u). If we insert ( f4(i| ) into ( f26j ) 
we get, at first order in n, 



6i(l + u) -bx{l-u) 



c I e 
la 



-c(v*-uv)/2 dy ( 41 ) 



It is easy to check that a solution of (41) compatible 



with the conditions &i(l) = and b\(u) = bi(—u) is 



h(u) 



+00 cosh ( A "/ 5 
dX 



cosh 







sinh 



(42) 

There are however many other solutions of (^l|), 
which can be obtained by adding to (^) an arbi- 
trary function F(u, c) even and periodic in u of pe- 
riod 2 and vanishing at u = 1. If we require that 
each term in the small c expansion of 61 (it) is poly- 
nomial in u (as justified in section 3T), we see that 
all the terms of the small c expansion of F(u, c) must 
be identically zero. This already shows that ( (42| ) has 
the same small c expansion ( |39| ) as what one would 
get by collecting all the terms proportional to n in 
the small c expansion of section |3.l| . 

If the solution ( 42 ) of (|l|) had a non-zero radius of 
convergence in c, it would be natural to choose this 
solution and set F(u,c) = 0. However it is easy to 
see that ( [i"2"| ) has a zero radius of convergence in c: by 
making the change of variable A 2 = 2^, it is easy to 
see that (^) is the Borel sum of a divergent series p5[. 

Apart from being the Borel sum of its expansion 
in powers of c, we did not find definitive reasons 
why (|4^) is the solution of (^lj) we should select. 
However, we can notice that for integer n, all the q a 
are real and B(u) defined by ( |24| ) is analytic in u and 
remains bounded as |Im u\ — > 00. The solution b\ (u) 
given by (|^j) is also analytic in u and grows as ln(w) 
as |Im it I — > 00. Adding any function F(u,c) peri- 
odic and analytic in u to ( 42 ) would produce a much 
faster growth. 



If we insert ([lO]) into (|2q), we have to solve at or- 
der n k 

b k (l + u)-b k (l-u) = ip k ( u )> ( 43 ) 
where (pk{u) is some function odd in u which 
can be calculated if we know the previous or- 
ders h(u), . . .,b k -i(u). 

fc — 1 „ u 

ip k {u) = cV / dve' c{v2 - uv)/2 h{l-v) x 

6*-<-i(l + u-t0. (44) 
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We see that the difficulty of selecting a solution of a 
difference equation appears at all orders in the expan- 
sion in powers of n, and we are now going to explain 
the procedure we have used to select one solution. 
If we write, as ipk(u) is an odd function of it, 

<p k (u) = 2^ dX sinh (^y^) a fc (A), (45) 

which is equivalent, by inverting when u is imaginary 
the Fourier transform in (|45|), to define a/c(A) by 

ak{x) = ^L dusin {^r) ^(.70' (46) 

then the solution for b^iu) we select is given by 
cosh ^ Atl v /5 ^ _ cosn 



/■+oo 

h(u) = / d\ • 
Jo 



sinh ( h& 



a k (X). 



(47) 



Indeed, bk(u) is an even function, vanishes at u — 1 
and one can check using (p5| ) that (|47|) solves ( 43j) . 

The integrals in (45|-|47|) are convergent JT7 and 
equations ( p4| , ff6| , |47 ) give an automatic way of cal- 
culating the bk(u) up to any desired order. 

This procedure is the direct generalisation of the 
choice (^2|) we did to solve (pd|). In fact, for fc = 
1, equations Q |46|) give (for A > 0) a x (A) = 
v /cexp(-AV2) and (^Jj) is identical to @. 

As for (|44 ), the solution (^) is not the only solu- 
tion of (|f3|). At any order fc, we could add an arbi- 
trary even periodic function F(u, c) of period 2, the 
expansion of which vanishes to all order in c. As for 
b\ (u), we did not find an unquestionable justification 
of our choice. One can notice nevertheless that (^) is 
the solution of ( f43| ) analytic in u and with the slowest 
growth when |Imu| — * 00. 



At order n , the procedure m4, Wq) gives 



02(A) = ce 



= r P ^l 2 



A 2 cosh 

tanh ' 



■(*) 



(48) 



d[i e 



-m 2 /2 



e -V/2 _ 2 



tanh 



(¥)J' 



with 62 (u) given by ( [f7j ) . Writing down 63 (u) or 03 (u) 
would take here about half a column. 

We can now give the first terms in the small n 
expansion of the energy. Using relation (129) , we find 



L 2 ( c c 2 \ 



(49) 



3/2 r + OO y2 



dX- 



-A 2 /2 



4 Jo tanhf^ 




'd^e-^/ 2 
'a tanh 

+ 0(n 4 ). 

By making the change of variable A 2 = 2i/, the terms 
of order n 2 and n 3 appear as Borel transforms of 
series in c with a finite radius of convergence. We 
conclude that these terms have both a zero radius of 
convergence in c. 

This small n expansion gives quickly very com- 
plicated expressions of bk(u). It turns out, as we 
shall see in the next section, that for large c, the 
expressions of the bk{u) get simpler and the en- 
ergy E(n, L, 7) can be calculated to all orders in pow- 
ers of n. 

4 Expansion in powers of n in 
the regime c — ► oo 

In the previous section, we have developed a proce- 
dure allowing to get the small n expansion of the en- 
ergy by solving the problem (|26|-|2lj|) . Here, we show 
how this procedure gets greatly simplified for large c. 

The expansion in powers of n of the previous sec- 
tion can be summarised as follows: if we use ( [I()| ) and 
we write 



a(A) = nai(X) + n 2 a2(A) + 



(50) 
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the bk (u) and dk (A) can be obtained by expanding in 
powers of n the following two equations 



B(u) = 1 



+00 cosh ( ^ ) - cosh 



It is apparent on fl5q ) and (|36|) that in the large c 
limit the function B(l + u/y/c) depends only on u 
and nVc, and a(A) depends only on A and riy/c. Let 
us introduce the constant K 



sinh(^) 



(this is a rewriting of (47)), and 



o(A) = 



nc 
7? 



/•+oo 




/ <iu sin 


:*)« 







o(A), 
(51) 

(52) 



K = 1 



Equation (55) becomes 



dAo(A). 



(57) 



B{l + - F )=K+ d\e Xu/2 a{\). (58) 



iu 



n 



V c 



(This is a rewriting of ( [44) , flq).) It will be conve- 
nient in the following to replace J5^ ) by its Fourier 
transform 



2 £Txsmh(^^)ai \ } 



(53) 



nc 



Jo 



(This is a rewriting of @ |||).) 

We are going to see how one can simplify (|5l|-|33|) 
when c is large. First we observe that for large c 
and u fixed of order 1, the expression bi(u) takes the 
scaling form 

r- f-oc 

6i(l + 4) - V~c / (e Au/2 ~ l)e- A2/2 dA. (54) 
V c Jo 

One can check from ([44], |^]) that this scaling 
form is present at any order in the small n expansion. 
Indeed, fl5l| ) becomes in the large c limit 



5(1 + -=) = 1 + / dA (e A "/ 2 - l)a(A), (55) 
Vc Jo 

and using (|5^), we find 

/>+oo / \ \ 

2 J dX sinh f a(A) = (56) 
nVcjdve-^-^B (i-J^\ B (l+^ 



In (|36|), if we write the integral from to u as the 
•difference between an integral from to +oo and an 
integral from u to +oo, and if we change the variable 
in the second integral to shift it to to +00, we obtain 



f Au 
2 dX sinh — o(A) = 
Jo 2 

/> + oo 

nVH / dve' v2 / 2 B{l-^) 
Jo V c 



(59) 



-uv/2 B(1 



U + V , 



If we replace B(l + (u — v)/y/c) and B(l— (u + v)/y/c) 
by their expression fl58| ) , we get after some rearrange- 
ments 



2 / dX sinh — a(A) 
Jo 2 



(60) 



/■+00 

^ / dve- v2 ' 2 B{\-^) 
Jo V c 

r+00 

/ a>a(^)e-^/ 2 2sinh 
Jo 



2ifsinhf^ 
V 2 



v + pi 



Taking the Fourier transform of this expression for 
imaginary u, we get for A > 



/•+00 

(A) =nV~c / dve- v2/2 B{l - -p) 
Jo Vc 



ii"5(A - «) + 



dAta(/x)e-' il,/2 5(A-t>-^) 



(61) 
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This last expression can be used to calculate -6(1 
u/yfc) using (||): 



K 



(62) 



V~c / dve~ v / 2 S(l-4=) 
■/o V c 



^a(M)e-' It,/2 e (l,+Al)u/2 



Finally, using (|58|), we recognise the relation 

/■+O0 

n^/c" / d ue -" 2 / 2 S(l - ^-)e OTl / 2 B(l + ■ 
■/o V c 



We see that, in the large c limit, ( pjj , p2\ ) reduce 
to this single equation ( |63| ) . We are now going to see 
that (|6^) can be solved to all order in the parame- 
ter n^/c. If we introduce the function (3{u) and the 
parameter e defined by 



(63) 



-)■ 



(3(u) 



2K^ 



(64) 



and 



(65) 



e = 2nK^/irc, 
then ( p3[ ) simply becomes 

1 /•+oo 

/?(u) = -^e" u / 4 + e / efo /3(u - v)p(-v). (66) 
2 V 7r Jo 



Using (g| _ 
energy i?(n, £, 7) in terms of (3{u 



E(n, L, 7) 



64 ) , we can express the ground state 

(3"(oy 



L 2 



n 3 c 2 



6 



nc- 



/3(0) 



(67) 



It is clear that relation ( |66| ) alone determines 
at least perturbatively in e. So, fro m (|67| ), we only 
need to extract /3(0) and /3"(0) from (|C>G[) . 

It is easy to do it for the first orders in e directly 
from equation (|36|). Moreover, we have found a way 
of calculating /3(0) and 0"(O), and hence the energy, 



to all orders in e. The details of the calculation are 
given in appendix |^. The final result can be written 

as 



1 



V 



v k=\ 



2 

E{n,L,i) = 



nc 
~12 



>c \ -» 



40F fc 5 / 2 
v fc=l 



(68) 
(69) 



We see that the energy is defined in an implicit 
way: expression ( pq) allows to calculate e as a func- 
tion of riy/c, and (169) gives the energy as a function 
of e. If we substitute c using (|l8|), we obtain the 
result announced in (|J). 

For small n, one can eliminate e from ( |68| ) and (|69|). 
We get 



^(n,L, 7 )-^-=^=(-_>^- (TO) 



^(2nV^) 2 



+ o ((nV^) z 




(2n v / CTF) 3 



5 Conclusion 

In this paper, we have calculated, using the replica 
method, the first cumulants ([[3], [l9|) of the free en- 
ergy of a directed polymer in a random medium ([!]) 
for a cylinder geometry. We used the integral equa- 
tion (|2|) of which together with conditions (§7j 
p8| ) allowed us to expand the moments of the 

partition function in powers of the strength c of the 
disorder or in powers of the number n of replica. All 
the coefficients of the small c expansion ( |38| ) are poly- 
nomial in n allowing to define the expansions for non- 
integer n. On the other hand, the coefficients of the 
expansion ( f49| ) in powers of n are complicated func- 
tions of c, with in general a zero radius of convergence 
at c = 0. As already mentioned in ]l7| ], we think 
that weak disorder expansions of the moments (Z n ) 
have generically a zero radius of convergence for non- 
integer n when the disorder is Gaussian; this is al- 
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ready the case for a single Ising spin in a Gaussian 
random field [jl7j. 

To obtain our small n expansion, we solved a dif- 
ference equation ( p6| ) which at each order in powers 
of n has several solutions. We selected the particular 
solution which has the slowest growth in the imagi- 
nary u direction and has the right small c expansion, 
but we could not exclude other solutions. A different 
approach, with a direct calculation of the first cu- 
mulants of the free energy, and not based on replica 
would therefore be very useful to test the validity of 
our expressions © which we have been able to de- 
rive only perturbatively to all orders in c. 

Although our expansion in powers of n becomes 
quickly very complicated, it simplifies when c is large 
and we could write in this limiting case, all the 
terms of the small n expansion (68, |69|). The ex- 



pression (H) we obtain of the energy E(n,L,j) (that 
is, through (||), the expression of (Z n )) is given ex- 
actly by the same scaling function as found for the 
ASEP. The present work therefore gives additional 
evidence that the scaling function G((3) given by (^], 
^) is characteristic of the long time behaviour of the 
KPZ equation in 1+1 dimensions on a ring and that 
the probability distribution of the free energy for a 
very long directed polymer on a ring should have an 
universal shape in the range where the fluctuations 
per unit length of the free energy are of order 1/L. 
Other universal distributions for the free energy of 
a directed polymer have been found recently for dif- 
ferent geometries p7 , 28, 2^, |30|| . Our present ap- 
proach, based on the Bethe ansatz is, at the moment, 
unable to recover these other distributions. One can 
try however to extend it to open boundary conditions 
(in this case too, the Bethe ansatz can be usedp^) 
instead of periodic boundary conditions and see how 
this change of boundary conditions affects the dis- 
tribution of In Z. Of course, it would be very nice 
to find a simpler approach which would somehow 
unify all these results and allow to relate all these 
universal distributions corresponding to the possible 
geometries, in the spirit of critical phenomena in two 
dimensions where conformal invariance^l] allows to 
connect the properties of different geometries. 

Technically, the approach followed in the present 
work is simply to try to find the q a solution of ( |l7| ) 



and to calculate the energy ( pf ) which is a symmetric 
function of the roots q a , in such a way that n becomes 
a continuous variable. One could do the same in all 
kinds of situations. For example, in appendix H we 
show how to define and calculate symmetric functions 
of the roots of Hermite polynomials when the degree 
of the polynomial becomes non-integer. 

Another interesting extension of the present work 
would be to consider more general correlations of the 
noise (^). The corresponding quantum problem be- 
comes then the general problem of quantum parti- 
cles interacting with an arbitrary pair potential. If 
the interactions are short ranged, one expects the 
universality class of the KPZ equation to hold, so 
one could try to repeat our expansion in powers of c 
for a general potential (without the use of the Bethe 
ansatz) simply by a standard perturbation theory in 
the strength of the potential. We believe that at any 
order in the strength of the potential, the ground 
state energy is polynomial in n allowing to define the 
perturbation expansion for non-integer n as we did 
here. If, with such an approach based on perturba- 
tion theory, one could recover the scaling function 
G of (|[ [?|), one could try to extend the approach 
to higher dimension as the relation between the di- 
rected polymer problem and the quantum hamilto- 
nian is valid in any dimension. 
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A Derivation of [27|, gS], gg) 



Let us first establish some useful properties of the 
numbers p(q a ) defined by (|25|). If the q a are the n 
roots of the polynomial P(X) 

P(X)=H(X-q a ), (A.l) 

it is easy to see that the p(q a ) defined in (|2^) satisfy 

(A.2) 



P(X + c) =1 + c y\P(9a) 



P(X) 



X 
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(The two sides have the same poles with the same 
residues and coincide at X — > oo.) Expanding the 
right hand side of (|A.2|) for large X, we get 



P(X 



P{X) 



= 1 



p(ga) 



gq , gq 

X X 2 



(A.3) 



On the other hand, using (|l6, AJ) and the symmetry 
{la} = {-9a} we have 



P(X) =X n + —E(n,L 7l )X n - 2 



+ o(x n - i ), 



so that 

PjX + c) 
P(X) 



nc 

= 1 + x 



X 2 



(A.4) 
(A.5) 



„3 / " 



+ 



3 )-aE(n,L, 7 )£ 2 /2 / 1 

I ^4 



X 3 



Comparing ( |A.3[ ) and (A.5), we get the relations 
n, 



^P(9a 
^9aP(fc) = c 2 (3) 



£(n,L, 7 )L 2 



(A.6) 
(A.7) 

(A.8) 



Moreover, by letting X = ±qg — c in (A. 2) we get for 
any qg root of P(X) 

1 = P(ga) = P(gq) 



9a - 9/3 + C 



9a + 9/3 + C 



(A.9) 



9a 9o 

Lastly using the symmetry {q a } = {— 9a} and the 
definition (|25|), the Bethe ansatz equations ( [Tt| ) re- 
duce to 



e 9a pHZa) -e- 9 «p(ga) =0. 



(A.10) 



From the definition ( p4| ) of -B(m) and the properties 
( [A!6HA.1C| ), it is straightforward to establish (§§-§!): 
the integral equation ( p6| ) is a direct consequence of 
(H) and jA~9|). Pro perties (|7|, |2|) follow from (||, 



consequence of <M4 



and (g4 A. 10) respectively. Lastly (|29|) is a 



A.6-A 



B The energy in the scaling 
regime 

In this appendix, we show how to calculate the energy 
from the integral equation (|66|). This equation is of 
the form 

r+oo 

(3{u)=H{u) + e / dvf3(u-v)(3{-v), (B.l) 



(B.2) 



where, in our case, H(u) is given by 

1 



H(u) 



2^ 



We are going to do our calculations for an arbitrary 
function H(u), even in u and decreasing fast enough 
(to make all the integrals converge) when \u\ — > 00. 

To find the energy, we see from (|67|), that we have 
to calculate from (B.l) the quantities /3(0) and /3"(0) 
as functions of e. We first show that (B.l) is equiva- 
lent to 



f3(u)=H(u) + e dvH(u-v)f3(v), (B.3) 
Jo 

as long as H(u) is even and decreases fast enough. 
Then, we will introduce a new function (3* (u) which 
is easy to calculate, and relate the derivatives of (3(u) 
and (5*(u) at u = 0. 



B.l Equivalence 
and (|B~3| ) 



between 



(ED 



The solution of ( p.3| ) can be written as 

I3{u) = p (u) + efa(u) + e 2 (3 2 {u) + . . . , (B.4) 
where 

A)(u) = H(u), (B.5) 
r+00 

0i(u)= / H(u-vi)H(vi)dui, 
Jo 
/• r+oo 

02{u) = J J H(u - vi)H(vi - v 2 )H{v 2 ) dvxdv 2 , 

/p+00 
■j H{u-v 1 )H{v l -v 2 )... 

■ ■ ■ H(v k -i - Vk)H(vk) dvi . . . dvk- 
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For a given k > 0, the integration range of /3k(u) can 
be divided into k parts: the region where v\ has the 
lowest value of all the the region where vi has 

the lowest value, . . . , the region where Vk has the 
lowest value. Let us consider, for some j such that 
1 < j < k, the region where Vj has the lowest value. 
All the other integrals then run from Vj to +00. If 
we translate those to integrals running from to +00 
by changing V{ into v t + Vj, we get: 



+00 
dvi 



+00 
dvi . 



+00 
dv 



3+1 



. dvj-i H(u — Vi — Vj) (B.6) 
H(vi -v 2 )...H(v j - l ) x 

...dv k H(-v j+ i) 

H(vj+i - v J+2 ) ■ ■ ■ H(vk + vj) 



Using the fact that H(u) = H(—u), we see that (B.6) 
is equal to 





By summing over j, we therefore have 



1% 



(tt) = / dv y^ y Pj^i(u - v)0 k _j(-v). 
■Jo J=1 



(B.7) 



(B. 



Finally, if we multiply by e k and if we sum over k all 
these terms (keeping apart the term for k = 0), we 



obtain equation (B.l 



Th e equatio ns (B.l) and (B__3) a re th us equivalent 
and (B.4, B.5) give the solution of (B.l) to any order 
in e. 

B.2 Calculation of the derivatives 

of p(u) 



If we look at the expression ( [B.5D of (3{u) in pow- 
ers of e, the calculation of (3(0) and /3"(0) looks sim- 
ple, especially when H(u) is given by ( p3.2[ ). How- 
ever, when we try to actually do the calculation, the 
expressions become quickly complicated with error- 
functions, primitives of error- functions, etc. It would 



be much easier if the integrals in (B.5) were running 



from —00 to +00 instead of to +00. This is why we 
introduce the even function 

(3* (u) = /3* (u) + eP* (u) + e 2 ff 2 («) + ..., (B.9) 

where, for k > 0, 



#(«) 



1 



k + 1 



H(u -vi) ... H(v k ) dvi .. . dv^ 

(B.10) 



and (3q(u) = H{u). One can see easily that 

P*{u) = ^ jdqe~ lqu ln(l - eH(q)), (B.ll) 
where we have defined 

/+00 
due lqu H{u). (B.12) 
-00 

The Wiener-Hopf technique (32| allows to re- 
late (3(u) and /3*(u). More specifically, we are going 
to show that for any X > 0, 



r+oc / r+00 

e Jdue- uX (3*{u) =ln(l + e Jdue- uX (3{u) 



(B.13) 



This relation allows to relate the derivatives of (3(u) 
and /3* (u) at u = 0: indeed, if X is large in flB.13 ), 
we get 



r e -- /3 ( U )=M + / 3 '(°) .mo) 



X 



x^ 



X 3 



(B.14) 



and a similar expression for (3*(u). Comparing both 
sides of (B.13) gives 



0(0) =/T(Q), (B.15) 
/?'(0) =|/?(0) 2 , 

0"(O) =/3*"(0) + ^/3(0) 3 . 

6 

(We have used the fact that /3*'(0) = because (3*(u) 
is an even function.) 
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In order to prove (B.13), the first thing to note 
is that, as H(u) decreases fast when u — > ±00, then 
also does (3(u). This allows to define the two "partial" 
Fourier transforms 



p+00 

$+(q)= due^(3(u), 
Jo 

$-(q) = fdue lqu p(u). 

J —OO 



(B.16) 
(B.17) 



It is easy to see that /?+(<?) is analytic in the upper 
half-plane (Imq > 0). Moreover, in this half-plane, 
$+(q) is bounded and vanishes when \q\ — > 00. Con- 
versely, (3- (q) is analytic, bounded and decreases to 
at infinity when Img < 0. 

The function (3(u) can be written in terms of P+(q) 
and (3-(q): 

0{u) = i (dqe~ tqu 0+{q) + £_(?))> (B.18) 



which allows to express the right-hand side of (B.13) 
when X is positive: 



+00 

du~- uX 



uX (3(u) 



27T J_* X + iq 



(B.19) 



1 

2^ 



X 



iq 



We calculate the two integrals in the right hand side 
of ( B.19 ) by the residue theorem. As (3+(q) is analytic 
and decreases at infinity in the upper half-plane, the 
first integral can be written when X > as a contour 
integral around the upper half-plane. The only con- 
tribution to the first integral comes, using Cauchy's 
theorem, from the pole q — iX. One can also check 
that the second integral vanishes (using a contour 
around the lower half-plane and the fact that (3-(q) 
has no pole). Therefore, (B.l£) gives 



+00 



due- uX fi{u) = (3+(iX). 



(B.20) 



Now, if we multiply ( |B.3| ) by exp(iqu) and if we 
integrate over u, we easily get for any real q 

(3+(q) + /L(g) = H(q) + eH{q)0+(q). (B.21) 



This relation between H(q), (3-(q) and (3+(q), to- 
gether with ( B.ll| ) gives 



0*(u) 



1 r + °° 

— dqe- t *[\n(l+e0+(q)) .13.22, 



ln(l-e/3_(<?) 



Using again that, in the upper half-plane, (3+(u) is 
analytic and vanishes at infinity, we see that, for a 
small enough e, the quantity ln(l + e(3+(q)) is also 
analytic and decreases to at infinity when Img > 
0. Similarly, ln(l — e(3-(q)) has the same properties 
for Imq < 0. This allows to calculate the left hand 
side of ( B.13 ) as we did for the right hand side. We 
find 

+00 1 

due- uX (3*{u) = - ln(l + e(3+(iX)). (B.23) 



Comparing (B.20) and (B.23) completes the proof 
of jB.13p. 



We can now give an expression of the energy. If we 
use the definition ( |B~j ) of H(u) in ( |B.ll| , [B 12| ), we 
find 



-« 2 /[4(fc+l)] 



This gives 

(3*(0) 

(3*"(0) 



I +°° £ fc 

2A§(fc + l) 3/2! 
^ +00 k 



(B.24) 

(B.25) 
(B.26) 



and, together with (B.15), these equations allows to 
give an expression of (3(0) and (3"(0). 
From ( p7| , p4| . p5| ), we see that 



e(3(0) = riy/c 
Then, using ( p3.15| ) we get 

e(3*(0) =ny/c, 



(B.27) 
(B.28) 



(3(0) 



nJc b 
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The energy is given by (|67|). We get: 



E(n,L,>y) = — 



^--cVS/r"(0) 



(B.29) 



And, finally, using relation ( B.25 , B.26 ), we ob- 
tain dum. 



C Hermite polynomials with a 
non- integer number of roots 

What we try to do in this whole paper is essentially to 
calculate ^ Q q 2 a (the energy) where {q a } is solution 
of (|l7l), in such a way that n appears as a continuous 
parameter. This allows us to obtain expressions of 
the energy for non-integer n. 

One can use the same procedure in other kinds of 
situations. A simple example which illustrates our 
calculations is the case of the zeroes of Hermite poly- 
nomials. 

The n-th Hermite polynomial H n (X) is the solu- 
tion polynomial in X with leading coefficient 1 of the 
differential equation p3 



1 



H'^X) - XH'JX) + nH n (X) = 0. 



(C.l) 



The polynomial H n (x) is of degree n and has the 
symmetry H n (X) = (—) n H n (—X). For example, we 
have Hi(X) = X 4 - 3A 2 + |. The n roots {h a } 
(1 < a < n) of H(X ) are real and distinct |34f| . 

By deriving ( |C.l| ) p times with respect to X, we 
see that, for all p, 

XH^\X) = \h^ +2) {X) + (n-p)H^(X). 

(C.2) 

This shows that the (n — p)-th Hermite polynomial 
is, up to a constant factor, equal to the p-th deriva- 
tive of H n (X). (This property will be used a lot in 
appendix 

Equation ( |C.l| ) can be used directly to calculate 
the first coefficients of H n (X) 



H n (X) = X n 



X 



n-2 



4 V4 



X 



n—4 



(C.3) 



Using ( JC.3] ), the symmetry of H{X) and the large X 
expansion 



H n (X) 



E 

p>0 



1 



(C.4) 



we can calculate the moments of the roots {h a } 
of H{X): 



1) 



n(n — 1) 



(2n-3), 



(C.5) 
(C.6) 



and so on. These moments are a priori defined only 
for integer n but as the expressions are polynomial in 
n, one can obviously extend their definition to non- 
integer n (similarly to what we do in the small c ex- 
pansion of B{u) in section (3.2)). 

To generate all the moments of the roots h a , it is 
convenient to consider the generating function 



Q(u) 



(C.7) 



which is reminiscent of the quantity P{u) 
defined in our quantum problem. (Us- 
ing ( p4| ) and ([34]) we can check that f3(u) cx 
exp(u^/c/2) J2 p(q a ) exp(g Q w/^/c).) 

The function Q(u) is hard to calculate for general 
n but we can expand it in powers of n. This can be 
done by considering 



H n (X) 



+ 00 



duQ(u)e- u 



x 



(C.8) 



which is defined only for X positive and large enough 
to make the integral converges. This function ^(X) 
is sol ution of a differential equation which follows 
from (C.l): 



1 



1 



-&(X) + -^(X f - X^(X) + n = 0. (C.9) 

To obtain an expansion in powers of n, we write 

*(A) =n*i(A) + n 2 * 2 (A) + ... (CIO) 
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Thus #ipO satisfies 



1 



- X* !(X) + 1 = 0. 



(C.ll) 



This differential equation can easily be solved, and 
the integration constant can be fixed using the re- 
quirement (JCTst) that, for large X, \Jf(X) ~ n/X 



+ 00 



*!(X) = due- uX - u /4 . (C.12) 



Then order of ( JC.9D gives 

i^(X)-X* 2 (X) + ^ 1 (X) 2 = 0, (C.13) 
the solution of which can be written as 

r+oc , + oo C Osh(-%) - 1 

W) = 2 due- uX ~ u ' A / dt ^ e" 4 . 

Jo Jo t 

(C.14) 

The procedure can be iterated to any order in n (of 
course expressions become more and more compli- 
cated). Using ( |C.8| ) and the expressions of 
and ^2(X) we can give an expression of Q(u): 



Q(u) = ne~ u2/A + 



2n 2 e' u2 / 4 / dt 



cosh 



(*) 



(C.15) 



-e -1 +0(n 3 ). 



Expanding this expression in powers of u , one calcu- 
late from this expression and from ( C.7 ) the terms 
linear and quadratic in n of all the moments of the 
h a - (The result s ag r ee fo r the second and the fourth 
moments with (|C.5|, |C.6|).) 



We noticed that for small n, the expression (C.15) 



corresponds to n roots h a distributed along the imag- 
inary axis with a Gaussian distribution. We do not 
know whether this is general and whether there ex- 
ists, for general non-integer n, a distribution of the 
roots h a in the complex p lane plan e which gives all 
moments calculated as in (|C.5| , |C.6| ) . 

It is interesting to notice the similarity between 
Q(u) and f3(u) defined in section [|. 



D The expansion in powers of c 
using Hermite polynomials 

In this appendix we show how to expand the solu- 
tion {q a } of ( |l7|) in powers of c for integer n. One 
can see from (jl7| ) that the roots q a scale for small c 
like y/c. It is thus convenient to rescale the polyno- 
mial P(X) defined in (|l9| ) and the q a in the following 
way: 



q a = r a y/c, 



(D.l) 



P{Xyfc) 



R(X). 



({r a } are thus the roots of R(X).) With these new 
variables, equation ( pc| ) becomes 

e r ^ c R{r a -V~c) + e- T ^ c R(r a + VH) = 0. (D.2) 

As the roots r a of R{X) are all distinct, this equation 
is equivalent to 



^R{X -yfc) + e- x ^ c R(X + yfc)= (D.3) 
2(coshXv/H + /(X))E(X), 



where f(X) is analytic (this follows from the fact 
that as R(X) is polynomial, f(X) defined by ( p.3| ) 
is obviously mcromorphic; moreover as the left hand 
side of ( |D.3| ) vanishes at all the roots of R(X), f(X) 
has no pole.) We are now going to solve (D.3) as a 
power series in c (i.e. find both f(X) and R(X) as 
power series in c). 



D.l Expansion of the polyno- 
mial R(X) 



We only have the single equation ( p3.3[ ) to obtain two 
quantities (R(X) and f(X)); however, using the fact 
that f{X) has no pole and R{X) is a polynomial, 
both quantities can be determined in a small c ex- 
pansion. Let us write 



R(X) = R (X) + cRiiX) + c 2 R 2 (X) 
f(X)= ch{X) + c 2 f 2 {X) + 



(D.4) 



where the fi(X) have no pole, Ro(X) is a polynomial 
of degree n (the term of highest degree in Rq{X) is 
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A") and all the Ri(X) (for i > 1) are polynomials 
of degree less than n. At first order in c, we find 



that (D.3) gives 



1 



R — XR' — fiRo 



(D.5) 



As fi(X) has no pole, it must be a polynomial. Be- 
cause Rp( X) is of degree n, we see by looking at both 
sides of ( |D.5| ) that, necessarily, A (A) = —n. We 



recognise then the differential equation (C.l) that de- 
fines Hermite polynomials. Therefore 



fi(X) 
Ro(x) 



-H(X). 



(D.6) 



We recover that way that the r a are the zeroes of 
the n-th Hermite polynomial when c is very small ]24fl. 



At next order in c, equation (D.3) gives 



^R" - XR[ + nR ± - f 2 H = 



(D.7) 



X*_ 
6 



H' 



X< 



H" 



X 
~G~ 



1 



(3) _ ± H W 

24 



As R\ and H are polynomials, (D/7) tells us that 
/2-ff is a polynomial too. We also know that /2(A) 
has no pole, thus it must be a polynomial. Ri(X) 
is of degree strictly less than n, so the expression 
R'l/2 — XRi + nRi is of degree strictly less than n. 
As H is of degree n, we recognise in ( D.7 ) an euclidian 
division of polynomials: —/2(A) is the quotient of the 
right hand side of equation (D.7) divided by H(X), 
and the terms involving R\{X) form the remainder 
of this division. This ensures that there is only one 



possible function /2(A) which verifies (D.7). 

In practise, to perform this euclidian division we 
can use the property ( |C.2] ) of the Hermite polynomi- 
als as many times as needed in the right hand side 
of ( D.7 ): for instance, we transform the term X 3 H' /6 
into nX 2 H/6 + X 2 H"/12. We cannot change X 2 H 
anymore, but we can apply ( |C.2| ) to the term X 2 H" . 
When no more transformation is possible, we are left 
with: 



X : 



H' 



X 2 
4 

r2 



H" 



X 



—X z — 
.6 



6 

n(n — 1) 



1 

2T 



(D. 



G 



H-—H". 
12 



The Euclidian division is then easy to perform 

n(n — 1) 



/ 2 (A) = -^A 2 + 



G 



(D.9) 



—R'{ — X R\ 
2 1 1 



1 



ni?i = H". 
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Using again (C.2), the differential equation on R\ can 
be solved; we find 

Ri(X) = -^H"(X). (D.10) 

As Ri(X) is simply a derivative of H(X), and as 
fi(X) is a known polynomial of X, we see that at 
the next order in c we will have to solve an equation 
of the form 



i?2 - XR 2 + nR 2 - f 3 H = ^A J # (fe) . (D.ll) 



1 



Using many times equation ( |C.2| ) the right hand side 
can be written in a "canonical form" : 



(D.12) 



which allows to write f$ as a polynomial in X and 
i?2 as a sum of derivatives of H(X). It is easy to 
see recursively that at any order c k in the expansion 
we can repeat this procedure to calculate fk(X) and 
i?fc_i(A). As a result we see that fu is a polyno- 
mial in X and that Rk-i can be written as a sum of 
derivatives of H{X). 

It is worth noting that at each order the variable 
n comes from the previous orders and from transfor- 
mations of the kind XH'(X) -> \H"{X) + nH(X). 
Because those are the two only mechanisms by which 
n appears, it is easy to see that at each order the 
coefficients of the sum of derivatives of H(X) that 
constitutes Rk—i(X) are all polynomials in n. 

A computer can easily do this tedious but straight- 
forward task to any desired order. Up to c 3 , we find: 



R=H- —H" - c 2 
24 



360 5760 



2520 3024 



31 



H" 



lln 
60480' 



967680 



ff (6) +0(c 4 ). (D.13) 
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D.2 Expansion of the roots r a of R(X) 

As seen in ( |D.13| ), the polynomial R(X) is to leading 
order in c given by H{X). It is thus natural to write 
the roots r a of R(X) as 



0(c 2 



(D.14) 



{\h a } are the roots of H). Inserting ( D.14| ) 
into (D.12), we find, at first order in c, 



x a H'(h c 



24 



H"{h a ) = 0. 



(D.15) 



Using the definition ( |C.l| ) of Hermite polynomials, we 
have H"{h a ) — 2h a H'(h a ). This gives in turn x a — 
jxh a . Repeating this procedure to any order in c, we 
generate terms of the form h J a H^ (h a ) which can be 
reduced to terms of the form h l a H'(h a ) by using ( p.2| ) 
as many times as necessary. It is then possible to di- 
vide the expression by H'(h a ) and we are left with 
an equation giving each new term in the expansion 
of r a as a polynomial in h a . Again, this can be pro- 
grammed, and we get, up to the order c 2 : 



9a _ , 

i— — '^a 



c 2 / / n 11 

12 a + C \ { 120 ~ 1440 



1 



(D.16) 



Using (D.l) and (iq), this leads to 



^(n,L, 7 ) = - C $> 2 a -^£^ 



(D.17) 



c 

360 



which coincides with 



when one uses the proper- 
ties (CJ3, |C.6[ ) of the roots h a of the Hermite poly- 
nomials. 
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